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Can microlensing fold caustics reveal a second stellar 
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ABSTRACT 

Dense high-precision photometry of microlensed stars during a fold-caustic passage can be 
used to reveal their brightness profiles from which the temperature of the stellar atmosphere 
as function of fractional radius can be derived. While the capabilities of current microlensing 
follow-up campaigns such as PLANET allowed for several precise measurements of linear 
limb-darkening coefficients, all attempts to reveal a second limb-darkening coefficient from 
such events have failed. It is shown that the residual signal of a second coefficient character- 
izing square-root limb darkening is ^ 25 times smaller which prevents a proper determination 
except for unlikely cases of very high caustic-peak-to-outside magnification ratios with no ad- 
equate event being observed so far or for source stars passing over a cusp singularity. Although 
the presence of limb darkening can be well established from the data, a reliable measurement 
of the index of an underlying power-law cannot be obtained. 
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1 INTRODUCTION 

Since Ischneider & Weld ^98^ have proposed to use caustics of 
gravitational lenses for measuring the brightness profile of closely- 
aligned background sources, dense high-precision photometry on 
microlensing events involving fold-caustic passages has provided 
meas urements of linear l imb-d a rkening coefficients for several 
stars jAIbrow et alj|200d 1200 it lAfonso et aljlioool: iFields et"ai] 
l2003h. 



2 INTENSITY PROFILE AND LIGHTCURVE 

With / denoting the average brightness and p the fractional radius 
of an observed source star, its brightness profile for a given filter 
reads 

iip)=iap), (1) 

where ^(p) is a dimensionless function fulfilling the normalization 



For stellar brightness profiles that involve linear, square, or 
square-root terms in cos i9, where i9 is the emergent angle, proper- 
ties of light curves during fold-caustic passages have been studied 
by^Rhie & Bennett ( 1999), while prospects and strategies for de- 
termining the linear limb-darkening coefficient have recently been 
discussed in detail bv Dominik (.2004.) , showing the ability of ob- 
taining precise measurements even with moderate use of the current 
capabilities. 

This letter will however show that the measurement of addi- 
tional coefficients, e.g. corresponding to a square-root law term, is 
not possible for typical events that involve a fold-caustic passage. 
Attempts of such measurements (Fields et al. 2003) have resulted 
in the meaningful measurement of a single linear combination of 
the involved limb-darkening coefficients only, while a similar re- 
sult ha s been obtained for an event where the source passes over 
a cusp ^Abe et al."'2003'), where however the exhibited differential 
magnification is modest due to the source size being relatively large 
co mpared to the cusp. In contrast, the cusp-passage event discussed 
bv lAlbrowetai]h999iJ) provides a better size ratio between source 
and cusp and to date presents the only measurement of both linear 
and square-root limb-darkening coefficients by microlensing. 



C(p)pdp= -. 



(2) 



Let us consider the source passing a fold caustic during the time- 
span 2 and let F[ denote its flux at the beginning of a caustic en- 
try or the end of a caustic exit at time tf. In the vicinity of the caus- 
tic, the observed flux can be approximated as fe.g. .Albrow et all 
rr999bl:lDominikl2003l) 



1/2 



(3) 



where upper (lower) signs correspond to caustic entries (exits), Fr 
denotes a characteristic rise flux, t is an arbitrarily chosen unit 
time, and the variation of images of the source that do not be- 
come critical as the caustic is approached is neglected. The caus- 
tic profile function Gf (77; £,) depends on the caustic passage phase 
77 = liz [t — t'^)/tj^ and is c haracteristic for t he adopted stellar 
brightness profile ^(p), where ^Hevrovsk^vl2003^ 
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and 



1/2 . / V- 1 



(4) 



(5) 



With K denoting the comple te elliptical integral of first kind, j{z) 
reads JOaudi & GoulJl999h 

for 2 < -1 

= { for -1<.<1 



3 POWER-LAW PROFILES 

A popular class of models for the stellar brightness profile ^(p) is 
a linear superposition of terms proportional to powers of cos i9 — 
•^1 — p^, where i9 denotes the angle between the normal to the 
stellar surface and the direction to the observer. With the term pro- 
portional to cos'' i9 reading 

^M(p) = (i+p/2)(i-p'r^', 



(7) 



the stellar brightness profile becomes 

k 



(8) 



with k coefficients ^5 ^{pi) ^ 1 corresponding to the power 
Pi > which moreover are required to fulfill the condition 



< 1. 



(9) 



According to Eq. J4}, the caustic profile function Gf (r;; 5) en- 
herits the linearity in F^p.} from ^(p), making it a superposition 

fe 

GK^;r{pi}...r{,,})= i + ^r{,,} [Gii^^y(^) - i] (lo) 

of the caustic profiles functions ^p.y{ri) = G( {rj; (^^p.y) which 
correspond to the stellar brightness profiles ^{p.}(p), where 



Glip}{ri) 



1 (1 + f)! 



max(l — ?7,1} 



^X + T] - 1 



dx. (11) 



max(l— ?7, — 1) 



4 MEASUREMENT OF SECOND COEFFICIENT 

The stellar brightness profiles corresponding to uniform brightness 
(p = 0), linear limb darkening (p = 1) and square-root limb dark- 
ening (p = 1/2) are shown in Fig.Q while the corresponding caus- 
tic profile functions are shown in Fig.|2| It can be seen that square- 
root limb darkening mediates between linear limb darkening and 
uniform brightness. 

Let us compare maximal square-root limb darkening corre- 
sponding to the caustic profile function Gf {i/2}('?) with its best 
approximation by a combination of uniform brightness and linear 
limb darkening 




0.2 0.4 0.6 0.8 1 

fractional stellar radius p 

Figure 1. Stellar brightnes.s profiles ^(p) corresponding to different lin- 
ear superpositions of power-law models. The solid lines refer to profiles 
^(p; r) = 1 — r + r (P) which are affine-linear in cos i?, where F = 
corresponds to uniform brightness, F = 1 to maximal linear limb darken- 
ing, and F = 0.541 to the best approximation in the sense of Eq. 1131 to 
maximal square-root limb darkening whose profile 5{i/2}(p) is shown as 
dashed line. 
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Figure 2. Caustic profile functions G'^{ri; that correspond to the stellar 
brightness profiles shown in Fig.lTI where fa: F) = (1— F) G*^jjj (r))+ 
F C* j-^ J (ry) are plotted as solid Unes for F = (uniform brightness), F = 
1 (maximal linear limb darkening), and F = 0.541 (best approximation 
to maximal square-root limb darkening), while the caustic profile function 
G* ^l^2}^''^^ ^™ maximal square-root limb darkening is shown as dashed 
line. 



Gnv;T) = {l-r)G*i^o)iv) + TG*n,y{v) (12) 

which depends on the choice of a norm on the function space. By 
requiring 



[Gf*(r?;r)-Gr,{i/2}M'd,7 



(13) 



to adopt a miminum, the approximation is in maximal accordance 
with the use of for sampled lightcurves and yields F = 0.541 
as optimal value. Figs.Qand|2|also show the stellar brightness pro- 
file ^(77; F) and the caustic profile function Gfij]; F), respectively, 
that corresponds to this value. Since the caustic probes a combi- 
nation of several fractional radii p for any caustic passage phase 
rj which is characterized by T{ri, p) according to Eq. ^4), signifi- 
cant differences between the brightness profiles are washed out in 



Can microlensing fold caustics reveal a second stellar limb-darkening coefficient? 3 



the corresponding caustic profile functions making the latter hardly 
distinguishable. 

Let us assume that the uncertainties a of the flux measure- 
ments F follow Poisson statistics, so that a = (F/Fl)^^^ , 
where Uf denotes the uncertainty for the flux Ff observed for 
the source being just outside the caustic. The ratio Z{rj) between 
signal-to-noise for flux deviations AF due to the different profiles 
and signal-to-noise for the flux Ff is then given by 



square-root vs. linear l.d. 



Z{i/2}iv) 



AF(r?)/a 



where 



"{1/2} 



(Jpeak/Jf*) — 1 
G'pcak,{l/2} 



(14) 



(15) 



with Fpcak denoting the flux at the caustic peak corresponding 
to the maximum of Gj ('?)' which is Gpeak,{i/2} ~ 1.431. 
For selected values of Fpcak/F[, -^{1/2} iv) is displayed in Fig.|3| 
where the linear limb-darkening coefficient has been optimized so 
that either the integrated squared deviation or the single largest ab- 
solute deviation is minimized. In general, the optimal choice of F 
depends on the sampling, and with only parts of the caustic pas- 
sage being sampled, the signal may further decrease, also affected 
by adopting best-fit values for t{ and tj^ that differ from the true 
values. While a uniform sampling is better represented by the for- 
mer choice for an optimal linear limb-darkening coefficient, the lat- 
ter choice better reflects the case where the beginning of a caustic 
entry or the end of the caustic exit is preferentially sampled being 
the region that provides the largest amount of information about the 
square-root limb-darkening coefficient. 

For all fold-caustic events observed so far, the caustic-peak- 
to-outside magnification ratio is Fpcak/^T J; 8, so that the mea- 
surement of r{i/2} with an absolute uncertainty of 0.1 from a 
single point requires a relative flux uncertainty outside the caus- 
tic CTf/Ff < 0.2% for typical events, where this limit roughly 
increases proportional to Fpcak/Jf*. 

For comparison, Fig.|4|shows the relative signal-to-noise ratio 
Z{r]) for the detection of linear limb darkening against uniform 
brightness, where 



Z{i}{v) = 0(1} 



Glioyiv) - G?,{i}(??) 
^"{i}G'f^{i}('7) + l 



and 



«{!} = 



(Fpeak/Ff*) 



G, 



peak, { 1 } 



(16) 



(17) 



with Gpcak,{i} ~ 1.470 being the maximum of Gf ^-^y (77). 

From a simulation of a caustic exit with Fpcak/Jf* ~ 12.5 
lasting 12 h and being sampled on average every 6 min with an ac- 
curacy (Jf/Ff — 1. 5 % in accordance wi t h the capa bilities of the 
PLANET campaign iPominik et all2002h . lDominiki t2004..) found 
that a linear limb-darkening coefficient F^ij = 0.5 can be mea- 
sured with a relative uncertainty of ~ 2.5 %. With the residual 
signal of an additional square-root limb-darkening coefficient be- 
ing ~25 times smaller, the absolute uncertainty in r^i/2} would 
be about 0.3 and not allow a meaningful measurement. In or- 
der to achieve an uncertainty of 0.1 for r{i/2} in an event with 
Fpca.k/Ff = 8, one would therefore require a relative uncertainty 
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Figure 3. Ratio .Z{i/2}('?) between signal-to-noise AF{T])/a for the 
detection of maximal square-root limb darkening against the best- 
approximating afflne-linear limb darkening and signal-to-noise / a{ for 
the flux measurement outside the caustic for selected caustic-peak-to- 
outside magnification ratios i^pcak/^f*- In order to obtain the best approx- 
imation, r has been chosen so that the integrated square of Z^i/2}{v) = 
{AF{ti) /a) / {Ff /a-{) is minimized for the curves shown in the upper 
panel, yielding T =0.542, 0.544, or 0.546 for Fp^ak/^f* =3.2, 8, or 20, 
respectively, whereas the absolute maximum of Z^i/2}{v) has been min- 
imized for the curves shown in the lower panel yielding F =0.567, 0.576, 



or 0.586 for the same values of F, 
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Figure 4. Ratio Z^^iy{'r]) between signal-to-noise AF(r])/cr for the de- 
tection of maximal linear limb darkening against uniform brightness and 
signal-to-noise F* / erf for the flux measurement outside the caustic for se- 
lected caustic-peak-to-outside magnification ratios Fpcak/Fj*. 



of the flux measurement outside the caustic of /Ff ~ 0.3 % and 
less at the peak. With systematic effects dominating the photomet- 
ric error bars rather than the photon noise at this level, this appears 
not to be achievable. In addition to large caustic-peak-to-outside 
magnification ratios Fpcak/^f*, large fluxes Ff at the caustic out- 
side also enhance the prospects for measuring limb-darkening co- 
efficients by allowing better photometric accuracies Of / Ff . 
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5 MIXED LINEAR AND SQUARE-ROOT LAWS 

With square-root limb darkening being fairly approximated by a su- 
perposition of linear limb-darkening and uniform brightness, pho- 
tometric observations during the passage of a source star over a 
fold caustic appear to be blind to the power index of an adopted 
limb-darkening law. 

Assuming the presence of both linear and square-root limb 
darkening, it is a linear combination of the corresponding limb- 
darkening coefficients that is probed in first instance. 

Let us consider that maximal square-root limb darkening 
^{1/2} (p) is best approximated by 5(p;7), where 7 depends on 
the data. The stellar brightness profile can then be written as 



C(p) = i + ^(f±|±(p)-i) 



(18) 



where 



+ 7 

i-ip) = [C{i/2}(P)-7C{1}(P)] (19) 

fulfill the normalization of Eq. J2j and analogous relations hold 
for Gf (77;^). The mixed profiles ^+{p) and ^-(p) have been 
chosen so that C+(p) combines best-matched equivalent profiles, 
whereas 5- (p) measures their difference. The corresponding limb- 
darkening coefficients are then given by 



r4 



1 + 7 
27 



[7r{i/2} + r{i}] 



f- = [7ro/2}-r{i}] 



(20) 



Therefore, the data will provide an accurate measurement of f + 
similar to F = r{i} if only linear limb darkening is considered 
(i.e. T^ipy = 0), where one finds the correspondence T+ = 
[(1 + 7)/(27)]r with 7 ~ 0.54 ...0.60. With the influence of 
^_ (p) on the lightcurve being ~ 25 times smaller than that of 
^+{p), the uncertainty in r_ corresponding to the second limb- 
darkening coefficient will exceed that in r+ by the same factor, 
and models that coincide in r+ are roughly equivalent. 

With measurements of r+ and r_, the coefficients corre- 
sponding to the linear and square-root term follow as 

7 f, 7 f, 



r{i} 

^{1/2} 



1 +7 ' 
1 

1 + 7 



r+ + 



1-7 
1 



darkening coefficient are ~ 25 times smaller, making its mea- 
surement with current microlensing follow-up campaigns such as 
PLANET impossible, unless the caustic-peak-to-exit magnification 
ratio becomes Fpcak/^f* ^ 40 for a (typical) photometric accuracy 
o-f / Ff = 1.5 % at the caustic exit where none of such events have 
been observed so far, or the source passes over a cusp singularity 
rather than a fold. 

If more than one limb-darkening base profile is assumed to 
contribute, only a single characteristic coefficient can be accu- 
rately measured which corresponds to a specific superposition of 
the corresponding base profiles as found for worked examples 
(Fields et al. 2003: Abe et al. 2003.). 
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6 SUMMARY AND FURTHER DISCUSSION 

If the stellar brightness profile ^ (p) is modelled by a linear superpo- 
sition of three base profile functions, as for the discussed examples 
of constant, linear, and square-root terms in cos i9, the normaliza- 
tion forces one of these profiles to mediate between the other two, 
so that this base profile can be approximated by a superposition of 
the others. With several fractional radii p being probed by the fold 
caustic for any passage phase according to T{ri, p), the residuals 
of this approximation are smaller for the caustic profile function 
G{{Tj; ^) than for the stellar brightness profile ^(p), which limits 
the power of photometric microlensing observations during fold- 
caustic passages for revealing the stellar brightness profile f (p). 

Compared to the determination of a linear limb-darkening co- 
efficient, the residual signals of an additional square-root limb- 



